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Abstract 

Wireless multi-hop ad hoc communication networks represent an infrastructure-less 
and self-organized generalization of todays wireless cellular networks. Connectivity 
within such a network is an important issue. Continuum percolation and technology- 
driven mutations thereof allow to address this issue in the static limit and to con- 
struct a simple distributed protocol, guaranteeing strong connectivity almost surely 
and independently of various typical uncorrelated and correlated random spatial 
patterns of participating ad hoc nodes. 
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1 Introduction 



Today's wireless communication mainly relies on cellular networks [1,2,3]. At 
first, the sending mobile device directly connects to its nearest base station. A 
backbone network then routes the communication packets to the cell, where 
the intended receiving mobile device is registered. Finally, the cell's base sta- 
tion transmits the passed-by message to the latter. As part of the centralized 
backbone infrastructure each base station acts as a router, possesses the net- 
work information, controls the single-hop communications within its cell and 
assigns different channels to its various mobile clients. The base stations need 
to be placed according to some optimized coverage layout. This requires an 
enormous planning effort ahead of operation and leads to a static infrastruc- 
ture, hard to change and adapt to new, revised needs. This costly inflexibility 
motivates a flexible and infrastructure- less peer-to-peer concept: selforganizing 
wireless mobile ad hoc networks [4,5,6,7]. 

In a wireless ad hoc network, a sending mobile device uses inbetween mobile 
devices to communicate with the intended receiver. Such a multi-hop connec- 
tion requires each mobile device to have additional router functionality. As 
a central control authority is missing, the participating devices need coordi- 
nation amongst themselves to ensure network connectivity, efficient discovery 
and execution of end-to-end routes and avoidance of data packet collisions on 
shared radio channels; of course, mobility of the devices also has an impact on 
the network performance, which has to be coped with. Contrary to these global 
network features, the selforganizing coordination rules, called protocols in the 
jargon of electrical engineers, have to be local. Due to its limited transmission 
range, a mobile device is able to communicate only with its current spatial 
neighbors. Hence, it can only extract information on its local surrounding. 
Since this is the only input into the coordination rule, the latter is by defini- 
tion local. Upon execution, it readjusts for example the device's transmission 
power to its new surrounding. 

In this Paper we focus on the important connectivity issue and ask: what 
is a good local coordination rule for transmission power management, which 
almost surely guarantees global connectivity for the whole network? We em- 
ploy a simple static model for ad hoc communication networks. This allows a 
connection to continuum percolation theory [8,9] based on random geometric 
graphs [10,11]. The spatially distributed ad hoc devices correspond to nodes, 
which are more or less locally connected by communication links. Two nodes 
establish a mutual link, only if the first node lies within the transmission range 
of the second and vice versa. For the case of constant, isotropic transmission 
ranges, a further mapping onto the classical picture of continuum percolation 
[12,13], stemming from the transport physics in continuous random media, is 
straightforward: whenever discs with radius equaling half of the transmission 
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range are placed around two nodes and overlap, the two nodes are linked. 

Continuum percolation allows to study, for example, the dependence of the 
probability for strong connectivity on the transmission range and to find a 
critical range, above which the ad hoc network graph is almost surely con- 
nected. The critical transmission range can be translated into a critical node 
neighborhood degree ngb cr \f Hence, a simple local coordination rule would be 
for each node to adjust its transmission power to yield a little above ngb cr \ t 
neighbors. As we will demonstrate, such a rule is not flexible enough to perform 
equally well in various different environments, like homogeneous vs. hetero- 
geneous random spatial arrangements of nodes or homogeneous vs. heteroge- 
neous propagation media. Fortunately, such a rule is able to give some guidance 
to develop local coordination rules with improved adaptation properties. In 
some respect, these new rules represent technology-driven mutations of the 
continuum percolation problem and demonstrate the usefulness to combine 
real-world needs in electrical engineering with modified concepts of statistical 
physics. 

In Section 2 a precise definition of random geometric graphs in the context 
of wireless ad hoc communication networks is given; various spatial point 
patterns, uncorrelated and correlated, are introduced. In Section 3 standard 
continuum percolation based on discs with constant and random radius is 
applied to ad hoc networks to obtain first rigorous statements on coordination 
rules for connectivity. A mutation of continuum percolation is presented in 
Section 4, which directly leads to a distributed, local coordination rule, flexible 
enough to cope with various environments. A conclusion and outlook is given 
in Section 5. 



2 Random geometric graph approach to ad hoc networks 

The topology of static ad hoc networks can be viewed as a random geometric 
graph. It consists of a spatial pattern of points, where each point is connected 
to some others by links. In Subsect. 2.1 the construction of three generic 
random point patterns is described. A simple propagation-receiver model is 
used in Subsect. 2.2 to establish links between points. Some selected geometric- 
graph features of interest are discussed in Subsect. 2.3. 

2. 1 Spatial point patterns 

Throughout this Paper three different generic random spatial point patterns 
are used: a homogeneous, a multifractal and a Manhattan point pattern. In 
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the following comes a short description of their construction. With no loss of 
generality a two-dimensional square area with side length L = 1 is used. 

In a random homogeneous point pattern each of the N points is given a random 
position (x,y) G [0, L\ x [0, L\. A typical realization is illustrated in Fig. la. 
By definition it does not show generic clustering. 

One way to construct simple clustered point patterns is to employ a binary 
multiplicative branching process. The nonuniform probability measure sup- 
ported on the unit square is constructed by iteration: at first the parent square 
is divided into four offspring squares with area 1/4. Two randomly chosen off- 
springs get a fraction (1 + /5)/4 of the parent probability mass fi — l, whereas 
the remaining two get a fraction (1 — f3)/A. In the next iteration step each off- 
spring square follows the same probabilistic branching rule and nonuniformly 
redistributes its probability mass onto its own four offsprings. After j iteration 
steps the probability mass /i = 1 has been nonuniformly subdivided onto 4 J 
subsquares with area 1/4P, where (j)2 J of these subsquares (0 < i < j) come 
with probability mass [(l+/3)/4] l [(l — /3)/4p'~\ One after the other each of the 
N points to be distributed is given an independent and uniform random num- 
ber between and 1, which, given some probability- mass- weighted ordering 
of the 4 J subsquares, corresponds to exactly one subsquare, onto the particle 
is deposited and randomly placed inside. One such realization of a point pat- 
tern is shown in Fig. lb. The hierarchical clustering of points is due to the 
hierarchical branching structure of the iteration process. - The probability 
measure constructed with such a multiplicative branching process is a mul- 
tifractal [14]. The construction of multifractal fields has some importance in 
such diverse fields as turbulence [15,16,17,18], finance [19,20], Internet traffic 
[21], high-energetic multiparticle dynamics [22] and deterministic chaos [23], 
just to name a few. 

As a third generic class of spatial point patterns a Manhattan street pattern is 
used. N x and N y streets are equidistantly placed parallel to the x- and y-axis, 
respectively. One after the other each of the N points is randomly placed onto 
one randomly chosen street. Fig. lc gives an illustration of one realization. 



2.2 Construction of communication links 

For a wireless communication network we define a link between two nodes i 
and j, only if they can communicate back and forth to each other. Let P = Pi 
denote the transmission power given to node i, then according to a simple 
propagation-receiver model, which does not account for shadowing and fast- 
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fading effects, node j is able to receive the signal once 
P/r a 

> snr . (f) 

noise 

r = Tij is the relative Euclidean distance between nodes % and j. The path- 
loss exponent a is assumed to be constant; for free-space propagation it is 
a — 2, but depending on specific in-/outdoor propagation it can vary typically 
between 1 < a < 6. For a successful signal transmission the received power 
Pi/rfj relative to a noise power needs to be larger than the minimum signal- 
to-noise ratio snr. Without any loss of generality the variables noise and 
snr can be set equal to one, implying a rescaling of the power P«. Condition 
(1) guarantees that node j is able to hear node i. This alone would define 
a directed link % — > j. Since bidirectionality increases the efficiency of the 
communication feedback, we focus on bidirectional links (i <-> j): also i needs 
to hear j, implying that also (1) has to be fulfilled with the substitution 
Pi — > Pj. We call a bidirectional link a communication link. 

The link construction can be given a simple geometric interpretation. With 
P = Pi the condition (1) translates into a maximum range R i = p}J a . Nodes 
that lie inside this circle with radius Ri around node % are able to hear this 
node. For a communication link to exist between nodes % and j, j has to lie 
inside i's circle with radius Ri and i has to lie inside j's circle with radius Rj. 
For the case that all Ri = R are identical, this link construction matches the 
standard link construction of continuous percolation [12]. 

Besides various spatial point patterns Fig. 1 also illustrates the communication 
links. Each node has been given identical power P/P noTm = 5, where the 
normalization P nor m = -Rnorm comes by setting 7iR^ orm = L 2 /N — 1/p equal to 
the reciprocal of the node density p. Note that no periodic boundary conditions 
have been used for this figure. 

2.3 Geometric- graph features of interest 

In Sections 3 and 4 various rules for assigning power values to the ad hoc 
nodes will be discussed. Each such rule together with a chosen generic class 
of spatial point patterns defines a specific ensemble of geometric graphs. For 
example, the graphs of Figs, la-c each represent one realization out of three 
different ensembles: (a) identical power for a random homogeneous spatial 
point pattern, (b) identical power for a multifractally clustered random spatial 
point pattern, and (c) identical power for a Manhattan random spatial point 
pattern. 

One question to ask for each ensemble of geometric graphs is for example: 
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how large is the average giant component? The giant component is defined 
as the magnitude of the largest connected cluster appearing in a graph real- 
ization; see again Fig. 1 for a visualization. A related and important question 
to ask for ad hoc communication networks is, what is the probability that all 
nodes are able to communicate to each other via multihop link routes? With 
other words, what is the probability that the giant component is equal to the 
number of nodes iV? As connectivity is for sure a very important issue for 
ad hoc communication networks, so is a generalization called /c-connectivity. 
A graph is called ^-connected if between every pair of nodes there exist at 
least k independent paths, which implies, that once k — 1 nodes are removed 
at random, the graph remains at least one-connected. A A;-connected ad hoc 
communication network is more flexible and robust to routing failure. Hence, 
another question: what is the probability for an ad hoc geometric graph to be 
/c-connected? 

A simple flooding algorithm is used to determine the giant component of an 
ad hoc graph realization: a random node is tagged in first place, then its 
neighbors are tagged, which then continue to tag their untagged neighbors, 
and so on, until the corresponding cluster is saturated. This procedure is re- 
peated for all untagged nodes, until all nodes of the graph are tagged. By 
definition, the largest found cluster is equal to the giant component. For an 
ad hoc geometric graph to be one-connected, i.e. strongly connected, the gi- 
ant component has to equal the total number of nodes. Another procedure 
to inquire one-connectivity uses the NxN Laplace matrix [24], which is the 
difference between the diagonal node degree matrix and the adjacency ma- 
trix. An element of the adjacency matrix is either one or zero, depending on 
whether a link does or does not exist between the two respective nodes; a 
diagonal element of the node degree matrix counts the link neighbors of the 
respective node. If the Laplace matrix possesses only one zero eigenvalue, then 
the graph is one-connected; the number of zero eigenvalues counts the num- 
ber of partitioned clusters. We have employed the Laplace matrix algorithm 
only as a supplement for small N geometric graphs. The probability for one- 
connectivity is estimated from a representative sample of ad hoc geometric 
graph realizations as the ratio between one-connected and all sample graphs. 

(/c>l)-connectivity is algorithmically very costly. In a nutshell, for each of 
the N(N — l)/2 pair of nodes belonging to one graph it has to be checked 
whether at least k independent paths exist. Although some faster, but ap- 
proximate algorithms can be found in the recent literature [25], we prefer to 
switch from costly k- to cheap pseudo-/c-connectivity. Given one-connectivity, 
the latter only requires each node to have at least k neighbors. In fact, a 
theorem exists [26] , which, translated into our language, guarantees for a geo- 
metric graph ensemble based on random homogeneous point patterns and the 
artificial constant-power link rule of Sect. 3.1 that in the large N limit the 
probabilities for k- and pseudo-/c-connectivity converge as they approach one. 
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3 Continuum percolation with artificial link rules 

3.1 Artificial link rule I: constant transmission power 

The simplest rule for power assignment is to allocate the same power value 
Pi = P to each node i. This is an artificial and unrealistic rule. It would re- 
quire all nodes either to be designed for only the same single-valued power 
operation or, given already network connectivity, to carry out fast synchro- 
nization; of course, also an outside provider could adjust all node powers to 
one single value, but this would give up the infrastructureless philosophy of 
ad hoc networks. Nevertheless, the constant-P rule is good to start with [27]. 

Fig. 2 shows the average relative giant component as a function of P obtained 
from a sample of 500 geometric graphs generated with random homogeneous 
spatial point patterns and the constant-P rule. The chosen rather small sample 
size produces already more than sufficient statistical convergence and keeps 
statistical error bars to such small values, that those will not be shown for this 
and the following figures. A percolation threshold behavior around P crit ~ 
4.5P norm is observed: for P <C P crit the average relative giant component is 
close to zero, whereas for P ^> P crit it is almost equal to one. The sharpness of 
the threshold depends on the number of participating nodes; with increasing 
N the transition becomes sharper. 

For the determination of the exact threshold position the limit N — > oo would 
be needed. To partially account for this, the relative giant component has also 
been determined from small and medium N simulations with periodic bound- 
ary conditions. Results are also shown in Fig. 2. The percolation transition 
becomes sharper and moves a little to the left once compared with the previous 
results, but still no full convergence for the employed increasing N values is 
obtained. Since for realistic ad hoc networks the limit N — > oo is out of reach, 
we do not pursue this matter further. For the rest of this Paper we focus on 
N = 1600, which accommodates best a hot-spot, i.e. big-crowd application of 
wireless ad hoc networks. 

The critical power P crit ps 4.5P norm , obtained by setting the path-loss exponent 
equal to a — 2, can be given a more illustrative interpretation: the factor 4.5 
reflects the average number of neighboring nodes, which is also denoted as 
the average node degree. Defining P crit = R 2 CT - lt along the lines of relation (1), 
we get (k) = nR 2 rit p ps 4.5irRf 10Tm p = 4.5. Picking one node, the probability 
to find k other nodes inside its circular disc with radius P cr j t is equal to 
p(jfe) = ( JV ~ 1 )g fc (l - q) N ~ 1 ' k « (A fc /A;!)e-\ where q = irR 2 crit /L 2 . For N large 
and q small, p(k) becomes a Poissonian with mean (k) — X — q(N — 1) ~ 
qN ~ 4.5. Note, that due to finite-size effects and the usage of no periodic 
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boundary conditions the actually sampled critical link degree is a little smaller 
than the asymptotic value k cr - lt = 4.53, which is stated for example in [10]. - 
This demonstrates that other than in terms of P cr i t it is also convenient to 
characterize the percolation phase transition in terms of R crit or /c crit . Once 
choosing other path-loss exponents a ^ 2, the former will change according 
to P ciit (a^2) = -R„7t 2 -Pcrit («=2) , whereas -R crit or /c crit remain as before. 

Next, spatial point patterns other than random homogeneous are discussed. 
Fig. 3 compares the relative average giant component as a function of trans- 
mission power obtained for random multifractal and Manhattan spatial point 
patterns with the random homogeneous case; consult again Fig. 1. 

Evidently for small P the relative giant component is larger for the multifrac- 
tal than for the homogeneous patterns, but for the convergence of the relative 
giant component towards one a much larger P is needed. Due to the hierarchi- 
cal clustering, subclusters of points are easily formed at small P since only a 
small transmission range is needed to connect the corresponding nodes. How- 
ever, in order to connect the various subclusters either directly to each other 
or via isolated nodes lying inbetween it needs a rather large power. Another 
consequence of the pronounced clustering is that the rather sharp percolation 
threshold observed for homogeneous point patterns blurs more the larger the 
splitting parameter f3 characterizing the multifractal point patterns is chosen. 

Also for the Manhattan point patterns the relative giant component increases 
faster for small transmission powers than for the homogeneous point patterns. 
The reason is that the average nearest-point distance is smaller for the ran- 
dom points confined to the one-dimensional Manhattan streets. The Manhat- 
tan threshold for the relative giant component is relatively sharp and comes 
at a transmission power value, which is somewhat smaller than for random 
homogeneous point patterns. - For consistency, it is instructive to map the 
continuous percolation based on random Manhattan point patterns onto the 
well-known square-lattice bond percolation [29]. Two nodes are assumed to lie 
on a one-dimensional straight line and to have a distance I = L/N x = L/N y 
corresponding to the distance of two successive Manhattan street crossings. 
They can only communicate with each other, if inbetween nodes come with 
successive distances x smaller than their transmission range R; otherwise the 
occurring void is too large to be bridged. For a one- dimensional Poissonian 
point pattern with density A = N/ (N x + N y ) the distance x of two consecu- 
tive nodes is exponentially distributed according to p(x) = Aexp (—Xx). This 
allows to estimate the probability for the occurrence of at least one too large 
void between the two picked nodes with distance /. The expression is 



p(void > R) = 

3=1 



rn 



(J-1)! 



(1-jRy- 1 l + -(l-jR) e 



-j\R 



(2) 
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where m = \l/R\; a derivation of this formula is given for example in Ref. 
[28]. Setting the path-loss exponent to a = 2 and according to (1) converting 
the threshold power P pa 4P norm into R we arrive at a value p(yoid>R) = 
0.47, which almost matches the critical bond probability Pbond= 0.50 of bond 
percolation on a square lattice [29]. 

As long as the average relative giant component is not exactly one, no direct 
knowledge on the probability for strong connectivity, i.e. one-connectivity, is 
possible. The (k—1) curve of Fig. 4 illustrates this quantity as a function of 
P obtained from a sample of 500 geometric graphs generated with random 
homogeneous spatial point patterns and the constant-P rule. It also reveals 
a threshold behavior, which sets in once the relative giant component has 
approached one. The relative factor of the respective threshold positions is 
about 2.3. The probabilities for pseudo-/c connectivity with k = 2 and 3 are 
also depicted in Fig. 4. Of course, their threshold is shifted to even larger P 
when compared to the pseudo-one threshold. - A part of Fig. 5 shows the 
probability for strong connectivity for random multifractal and Manhattan 
point patterns. Whereas the Manhattan curve is close to the homogeneous 
curve, the onset for a nonvanishing probability in case of the multifractal 
point patterns is shifted to extremely large values of the transmission power. 
This is due to the inhomogeneous spatial clustering and demonstrates that 
the constant-P rule is not efficient in such an environment. 

3.2 Artificial link rule II: iid transmission power 

So far all nodes were assigned the same transmission power P. This constant- 
P rule is not able to counterbalance spatially sparse regions, which remain 
disconnected to the giant component. Some long-range links are called for to 
establish connections between otherwise separated subclusters. Some ad hoc 
nodes then have to send with a transmission power larger than average. The 
simplest heterogeneous rule in this context is the iid-P link rule [30]. It treats 
the transmission power values assigned to the ad hoc nodes as independently 
and identically distributed (iid) random variables. Independently from the 
other nodes each node chooses its power according to the same probability 
distribution p(P) with mean (P). Translating P into a transmission range R, 
this rule places a disc with random radius R around each node. 

As a flexible representative a bimodal distribution 

M = p^w/( p ~ - *)) + ~ W + ft» ( 3 ) 

is chosen. It comes with two parameters j3i and #2 determining the variance 
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and skewness of this distribution, but leaving the mean (P) untouched. The 
previously used constant-P rule is reproduced once f3\ = or f3 2 = 0. Fig. 6 
illustrates simulation results obtained with random homogeneous point pat- 
terns. The relative giant component as a function of the average transmission 
power (P) is shown for some combinations of (5\ and /?2- All settings result in 
a shift of the percolation threshold to larger (P) values when compared to the 
outcome of the constant-P rule. Without showing we remark that similar re- 
sults are obtained for other point patterns, i.e. multifractal or Manhattan, and 
other power distributions, e.g. of lognormal type. The randomness and spatial 
decorrelation of the long-range links introduced by the iid-P link rule does 
not allow for a shift of the percolation threshold towards smaller (P) values. 
In the next Section a much more elegant and spatially correlated approach is 
found to achieve this goal. 



4 Continuum percolation with a local link rule 

Since a central control authority does not exist in pure ad hoc networks, there 
is no external provider to assign specific power values to the ad hoc nodes ac- 
cording to its global rule. It is the ad hoc network by itself which has to decide 
which power values are assigned to its participating nodes. Due to the finite 
communication range of each node, see again Eq. (1), these coordination rules 
have to be local. In the jargon of electrical engineers local means distributed. 
Such a local rule will be presented and discussed in this Section. 

4-1 Local link rule: minimum degree 

By exchanging so-called hello and hello-reply messages each ad hoc node is able 
to access direct information only from its immediate neighbors, defined by its 
links. The simplest local observable for a node is the number of its links, which 
is equal to the number of its one-hop neighbors. Based on this observable alone, 
a simple strategy for a node would be to decrease/increase its transmission 
power once it has more/less than enough neighbors [31]. Consequently the 
target node degree would be confined between a lower and upper bound ngb 
and ngb; for simplicity of the following arguments we set ngb = ngb = ngb. 
A value of the latter has to be chosen such that for example almost all nodes 
are part of one connected network and reflects the only external input to 
this otherwise local link rule. For example, the results of the previous Section 
obtained with random homogeneous point patterns suggest, that in order to 
guarantee a probability almost equal to one for one-connectivity the target 
value should be of the order 0{ngb) 15—20; consult again Fig. 4 and note, 
that according to the argument given in Sect. 3.1 in conjunction with random 
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homogeneously distributed point patterns and a path- loss exponent a — 2, the 
relative transmission power Pj P norm can be interpreted as the average number 
of node neighbors. 

This simple ngb local link rule has at least two drawbacks. The target range 
0(ngb) pa 15—20 might be sufficient for a randomly homogeneously patterned 
world of points, but for some other underlying point patterns it would only 
yield a probability for one-connectivity a little above zero. For the random mul- 
tifractal case, depicted in Fig. 5, the relative power has to be P/P norm > 40 
for this probability to come close to one, which corresponds to the numerically 
determined value 0(ngb) > 52. Besides this sensitivity on the specific nature 
of the point patterns, also the target values 0(ngb) ~ 15—20 and above are 
technologically unwanted, since they lead to too much blocking for the shared 
medium access control avoiding data-packet collision. The second drawback 
leads to frustration due to the specific nature of the bidirectional link defini- 
tion. A cluster of spatially close-by nodes might saturate, meaning that each 
of the nodes has ngb neighbors. A new node, located not so far away from 
this cluster, wants to connect to some of its nodes. In fact, the nodes of the 
cluster are able to hear the lonely node due to its large transmission range, but 
since they are already saturated they do not increase their power to bridge the 
necessary spatial distance. The links are only one-directed, but not bidirected. 
Due to the missing feedback, the lonely node further increases its power, even- 
tually up to its upper limit P max , remaining in a frustrated state of having too 
few neighbors and unintentionally interfering the others' communication. 

In order to avoid these drawbacks we present a modified local link rule. Upon 
setting up the communication links to the other nodes, a node attaches to 
its hello message information about its current link neighborhood list and its 
current transmission power. Starting with P m i n , the node increases its trans- 
mission power by a small amount once it has not reached a minimum link 
degree ngb min . Whenever another node, which so far does not belong to the 
neighborhood list, hears the hello message of the original node for the first 
time, it realizes that the latter has too few neighbors, either sets its power 
equal to the transmission power of the hello-sending node or leaves it as be- 
fore, whichever is larger, and answers the hello message. Now the original and 
new node are able to communicate back and forth and have established a new 
link. The original node adds one new node to its neighborhood list. Only once 
the required minimum link degree is reached, the original node stops increas- 
ing its power for its hello transmissions. At the end each node has at least 
ngb min neighbors. Some have more because they have been forced to answer 
nodes too low in ngb; their transmission power is larger than necessary to 
obtain only ngb m i n neighbors for themselves. 

Fig. 7 illustrates the algorithmic implementation of the local minimum-link- 
degree rule in more detail. Initially, all nodes come with a minimum transmis- 
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sion power Pj = P m j n and an empty neighborhood list Mi = 0. All of them 
start in the receive mode. Then, at random, one of the nodes switches into the 
discovery mode. By subsequently sending hello messages and receiving hello 
replies, the picked node increases its power until it has discovered enough 
neighbors. Then the node returns into the receive mode. For simplicity we 
assume that only one node at a time is in the discovery mode; furthermore, 
we assume the maximum transmission power P max to be sufficiently large, so 
that each node is able to discover at least ng& min neighbors. Another node, 
which enters the discovery mode at a later time, performs the same opera- 
tions. If during its previous receive-mode period this node had already been 
sending hello replies to then discovering nodes, it has to execute one more op- 
eration before returning into the receive mode: it compares its final discovery 
transmission power Pf sc with the maximum power Pj it has been asked to 
transmit hello replies; if the former is smaller than the latter, the node up- 
dates its neighborhood list by sending a further hello message with the power 
max(Pf lsc , Pi) and receiving additional hello replies. - In the receive mode a 
node listens to incoming hello messages. Upon receipt of such a message, the 
node first checks whether it already belongs to the incoming neighborhood list. 
If yes, the requesting node has already asked before with a smaller discovery 
power and there is no need for the receiving node to react. Otherwise, it up- 
dates its transmission power to max(P,, Pj), but only if the magnitude of the 
incoming neighborhood list is smaller than the required minimum link degree. 
Then it sends back a hello reply. If the node in the receive mode has already 
executed its discovery mode before, an additional operation is needed to up- 
date its neighborhood list: it sends again a full hello message and listens to the 
newly triggered hello replies. Upon doing so, the node checks on the status of 
the initial hello reply, whether node j was able to receive it or not, and, if yes, 
makes sure that both nodes are registered in their mutual neighborhood lists. 
Also some hidden links are identified by this additional procedure; a hidden 
link occurs when two nodes, having already executed their discovery mode, 
are forced to increase their respective transmission powers by an independent 
third and fourth party in such a way, that they are then able to communicate 
directly to each other, but do not yet have this knowledge. 

The implementation of the local minimum-link-degree rule written down in 
Fig. 7a represents the ad hoc node's view, which the latter use to selforga- 
nize into a network. The transmission power of each node is locally chosen to 
adapt to the spatial surrounding and is not globally assigned from outside. - 
As simulators we employ an equivalent, but much simpler from-outside imple- 
mentation. It is depicted in Fig. 7b. For each node % its relative distances 
to the other nodes j of the point pattern are sorted in increasing order. The 
first ngb m i n nodes of this list make up the discovery list Mf lsc . Invoking (1), 
the discovery power Pj = Pf lsc = (Vij(n S f> min )) a is fixed by the relative distance 
belonging to the n<?& m i n 'th node of the sorted list, but needs to be updated 
according to p = max(P, (r^)") once the node i falls into the discovery neigh- 
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borhood of the other nodes j. After completion of the heterogeneous power 
assignment, the links are easily constructed along the lines of the definition 
given in Sect. 2.2, which then fix the final neighborhood lists A^. 



4-2 Giant component and one- connectivity 



Figs. 5 and 8 illustrate the simulation results obtained with the local minimum- 
link-degree rule and compares them to the respective outcomes of the artificial 
constant-P rule. For geometric graphs based on random homogeneous spatial 
point patterns the threshold of the average giant component is reduced by 
about a factor of 1.45; for the probability for one-connectivity the threshold 
reduction factor around 2.05 is even slightly larger. Note that for ngb min = 3 
the relative giant component is already very close to one; once ngb m i n > 6 the 
probability for one-connectivity becomes one almost surely. 

For geometric graphs based on random multifractal spatial point patterns the 
local minimum-link-degree rule beats the artificial constant P rule even more 
impressively. Whereas for the artificial rule the giant component threshold is 
very blurred, the local rule transforms it into a sharp threshold, which almost 
exactly coincides with the respective threshold obtained for the previously dis- 
cussed random homogeneous point patterns. Due to the strong spatial cluster- 
ing the artificial constant-P rule leads to a highly suppressed one- connectivity 
yield; for the parameters used for Fig. 5 the onset for nonvanishing proba- 
bility is already at a rather large transmission power around P/P nor m ~ 13, 
but it needs an even much larger P for this probability to come close to one. 
The local rule, on the other side, acts like a wondrous outperformer, compen- 
sating again the strong spatial clustering introduced by the inhomogeneous 
spatial point patterns and pushing the one-connectivity threshold down to 
(P)/P nor m ~ 6, nearly matching the respective homogeneous-point-pattern 
threshold; with ngb min m 7, which is equivalent to (P)/P n0 rm ~ 10, the prob- 
ability for one-connectivity is practically one. 

In case of the geometric graphs based on random Manhattan spatial point pat- 
terns and in comparison with the artificial constant-P rule, the local minimum- 
link-degree rule reduces the threshold of the average giant component as well 
as of the probability for one-connectivity by about a factor of 1.7. The magni- 
tude of reduction is comparable to the values stated for random homogeneous 
point patterns. For the one-connectivity probability to become almost one, a 
value of at least ngb min = 10 is needed for the model parameters stated in Fig. 
5. 
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4-3 Node degree and power transmission distribution 

Next, we illustrate the node degree and transmission power distributions 
resulting from the minimum-link-degree rule and compare them with their 
constant-P counterparts. At first random homogeneous spatial point patterns 
are considered and the minimum link-degree is chosen to be ngb m i n = 6 in 
order to guarantee one- connectivity almost surely; consult again Fig. 5. The 
node degree distribution p(k), which reflects the probability for a node to have 
k neighbors, is shown in Fig. 9(al). By construction, p(k) = for k < ng& min . 
It comes with a peak at k — ngb m i n and falls off sharply for k > ngb m i n . 
Upon changing from discrete to continuous k > ngb m i n , the distribution can 
be nicely fitted as the superposition p(k) = a5(k — ngb ni \ n ) + bN{k; fi, a) of 
a (^-function, placed at the minimum link-degree, and a normalized Gaussian 
with shift /i and width a. The parameters used for the fit in Fig. 9(al) are 
a = 0.14, b = 1.49, fi = 5.93 and a = 2.39. - For comparison we also show 
the much broader node degree distribution resulting from the constant-P rule. 
The setting P/P nor m = 20 is necessary to guarantee one-connectivity almost 
surely; consult again Fig. 5. This distribution can be nicely fitted with a nor- 
malized Gaussian; parameters used in Fig. 9(al) are /i = 18.98 and a = 4.92. 
The small deviations to the expected Poissonian p(k) = (\ k /k\)e~ x with mean 
A = P/P n orm (see Section 3.1), which is also illustrated in this Figure, are due 
to finite-size effects, that on average nodes close to the boundary experience 
a degree less than A. 

The power transmission distribution in case of the constant-P rule is simply 
a 5-function. It is indicated as an arrow within Fig. 9(a2). Its position nat- 
urally represents an approximate upper bound for power transmission values 
obtained from the minimum-link-degree rule. For the distribution resulting 
from the minimum-link-degree rule an analytic estimate can be given. For 
simplicity the path-loss exponent is set to a = 2. As discussed already in Sect. 
3.1, the degree distribution of a node, not too close to a boundary, is given 
by the Poissonian p(k) = (qN) k e~ qN /k\, with q = (P/P norm )(7rP norm /L 2 ) = 
(P/Pnorm) /N representing the area covered by the node's transmission radius 
relative to the total domain area L 2 — 1. Here, the transmission power P is 
kept fixed and the degree k is the discrete random variable. Another look on 
Fig. 9(al) reveals, that the very narrow degree distribution resulting from the 
minimum-link-degree rule can be crudely approximated as p{k) pa 5{k — k ) 
with k = (k) = 7.36. Hence, the node degree is now kept fixed within the 
above Poissonian and the transmission power is considered as the continuous 
random variable. This leads to the transmission power distribution 

p(P/Pnorm) ~ (P/Pnorm)" eX P (-P/P norm ) . (4) 

It corresponds to a Gamma distribution p(x; a, b) = x a - 1 e- x / h /(b a T(a)) with 
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x = P/Pnorm and a = ko + 1 ~ 8.36, 6=1. The actual best fit to the sampled 
distribution, shown in Fig. 9(a2), yields the parameters a = 7.74, 6 = 1.01 and 
more or less confirms the given estimate. A fit with a lognormal distribution 
p(x) = exp{ — (lnx — /i) 2 /2<T 2 }/(v / 27T(Ta;) is also illustrated; parameter values 
are fi = 2.02 and a = 0.37. 

Node degree and transmission power distributions in connection with random 
multifractal spatial point patterns are illustrated in Figs. 9(bl) and (b2). The 
minimum link-degree rule produces a node degree distribution, which as for 
the homogeneous example can be nicely approximated by a superposition of 
a 5-function and a normalized Gaussian for k > ngb ni { n = 7. The parameters 
used for the fit in Fig. 9(bl) are a = 0.13, 6 = 1.36, // = 7.60 and a = 3.27. 
For comparison, the node degree distribution obtained from the constant-P 
rule is also illustrated. In order to guarantee one- connectivity almost surely, 
a rather large P/Pnorm = 50 has to be picked. This leads to a very broad 
distribution, which comes with an average degree (k) = 62.7 and which can 
be represented as a superposition p(k) = aiN(k; fj,i, G\) + a 2 N(k; /j, 2 , cr 2 ) + 
asN(k; /X3, 03) of three normalized Gaussians with parameters a\ = 0.04, \i\ = 
17.8, ai = 5.5, a 2 = 0.48, /i 2 = 40.4, a 2 = 13.7, a 3 = 0.48, // 3 = 89.5, 
0-3 = 22.0. The occurence of the double-hump structure is specific to the chosen 
parameters N = 1600, (3 = 0.4 and P/P norm = 50, where the transmission 
range area associated to the transmission power smoothes out any substructure 
generated beyond the iteration step j ~ 2—3 of the multifractal point pattern 
construction. - The sampled transmission power distribution resulting from 
the minimum link-degree rule is exemplified in Fig. 9(b2). A best fit to a 
Gamma and a lognormal distribution reveals better agreement with the latter. 
Best parameters are a = 2.11, 6 = 3.94, and ji = 1.97, a = 0.81, respectively. 

For completeness, node degree and transmission power distributions in con- 
nection with random Manhattan spatial point patterns are illustrated in Figs. 
9(cl) and (c2). The choice P/P n0 rm = 25 guarantees one connectivity al- 
most surely in context with the constant-P rule and leads to a rather broad 
Gaussian-like node degree distribution. The corresponding curve in Fig. 9(cl) 
comes with mean \x = 28.3 and width a = 6.8. Note, that the transmission ra- 
dius corresponding to P/P n0 rm = 25 equals about half the length 1/N X = 1/N y 
of a Manhattan block, so that the influence of the Manhattan structure still 
impacts the node degree distribution and has not been washed out to coin- 
cide with a node degree distribution of a random homogeneous point pattern 
with /i Rs 25. This difference also holds for the node degree distribution re- 
sulting from the minimum link-degree rule. The peak at k — ngb min = 10 
is very pronounced. The overall mean is (k) = 13.53 and the distribution 
p(k > ngb min ) w a8(k — ngb m[n ) + biN(k; /xi, o"i) + b 2 N(k; fi 2 , a 2 ) can be ap- 
proximated as a superposition of a ^-function and two normalized Gaussians. 
Parameters used in Fig. 9(cl) are a = 0.35, 61 = 0.63, fii = 7.78, <j\ = 2.31, 
and 6 2 = 0.56, /i 2 = 15.90, a 2 = 4.60. - The sampled transmission power dis- 
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tribution resulting from the minimum link-degree rule is exemplified in Fig. 
9(c2). A best fit to a Gamma and a lognormal distribution reveals a slightly 
better agreement with the latter, but deviations from both are clearly visible. 
Best parameters are a = 4.00, b = 2.57, and ji = 2.27, a = 0.54, respectively. 



5 Conclusion 

For wireless mobile ad hoc communication networks connectivity represents 
an important issue. In a selforganizing manner, the participating ad hoc nodes 
have to tune their transmission powers to establish direct one-hop communi- 
cation links to their spatial neighbors and to be able to reach all others via 
multihop routes. In the static limit, two-dimensional continuum percolation 
based on discs with constant or random radius can be mapped onto artificial 
link coordination rules with constant or random transmission power. Basically 
the transmission powers have to be chosen above a percolation threshold in 
order to guarantee strong network connectivity almost surely. However, the 
percolation threshold does show a sensitive dependence on the specific spatial 
patterning of the ad hoc nodes. Different classes of uncorrelated and corre- 
lated random point patterns, like homogeneous, multifractal or Manhattan- 
like distributions, make a big difference. The artificial link rules are not flexi- 
ble enough to adapt to local spatial inhomogeneities. A local generalization of 
these rules leads to the minimum-link-degree rule. It can be viewed as a step 
towards a selforganizing mutation of continuum percolation. It requires each 
ad hoc node to be connected to a minimum number of closest neighbors, which 
sets a lower bound on the node's transmission power. This node, on the other 
hand, might be forced to increase its transmission power further, in order to 
establish an additional communication link to another node, which belongs to 
the greater vicinity and has not yet reached the required minimum number 
of closest neighbors. This distributed rule is able to counterbalance local spa- 
tial inhomogeneities occurring in the random point patterns. As a function of 
average transmission power the percolation thresholds associated to the var- 
ious considered classes of random point patterns almost collapse onto each 
other and are tremendously reduced, when compared to the outcomes with 
the artificial rules. 

Compared to other local link rules yielding strongly connected networks, the 
presented minimum-link-degree rule is simple. For example, the Rodoplu & 
Meng algorithm [32] relies on GPS-based position knowledge of the ad hoc 
nodes to construct a link enclosure for each node; the construction also re- 
quires knowledge about the assumed uniform path- loss exponent, which can 
be seen as another drawback of this algorithm. Another proposal [33] requires 
ad hoc nodes to come with directional antennas; strong network connectiv- 
ity is then guaranteed once each node has neighbors in each angular sector. 
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Rules like these require more technological equipment for the ad hoc nodes. 
In terms of the simplicity principle, the presented minimum-link-degree rule 
appears to be more attractive. It also has the advantage that it will work in a 
heterogeneous propagation medium, where the path-loss exponent is a func- 
tion of position, distance and direction. A generalization of the distributed 
minimum-link-degree rule, which so far only has been developed for the static 
limit, towards mobile ad hoc networks is straightforward. These last two issues 
will be discussed in more detail in future work. 

Finding efficient distributed coordination rules, i.e. protocols, for connectivity 
is certainly one important issue for ad hoc communication networks, but there 
are definitely also others: power consumption, efficient routing discovery and 
execution, medium access control, interference, quality of service, and end- 
to-end throughput. The construction of an optimized protocol for one part 
alone, needs not be the overall best for all of the partially conflicting entities 
considered together. Clearly, the design of a rather complex, but still simple 
overall protocol is called for. It is possible that such an optimized protocol 
might lead to small-world or scale-free geometric-graph topologies, which are 
already observed and discussed for Internet and biochemical communication 
networks [34,35]. 
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Fig. 1. Geometric graphs for random patterns of N = 200 points confined to a box: 
homogeneous (top), multifractal (middle) and Manhattan (bottom). Parameters of 
the point patterns are j = 5, (3 = 0.4 (multifractal) and N x = N y = 7 (Manhattan). 
Each point is given the same transmission power P/P norm = 5, corresponding to a 
link range R/L = 0.089 when using the path-loss exponent a = 2. Points connected 
by solid links belong to the giant-component cluster. 
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Fig. 2. Average relative giant component as a function of transmission power. 
A sample of 500 geometric graphs generated with random homogeneous spatial 
point patterns and the constant-P rule has been used. The path-loss exponent 
of Eq. (1) has been set to a = 2. Different curves correspond to different num- 
ber of nodes: N = 100 (dash-dotted), 400 (dashed), 1600 (solid); curves marked 
without/with open circles correspond to exclusion/inclusion of periodic boundary 
conditions (pbc). 
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Fig. 3. Average relative giant component as a function of transmission power upon 
using the constant-P rule. The number of nodes has been fixed to TV = 1600 and 
the path-loss exponent has been set to a = 2. Different curves correspond to dif- 
ferent random spatial point patterns: homogeneous (solid), multifractal (dashed) 
with parameters (3 = 0.4 and j = 5, and Manhattan (dash-dotted) with parameters 
N x = N y = 7. A sample of 500 geometric graphs has been used for each case. 
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Fig. 4. Probability of pseudo /c-connectivity as a function of transmission power 
upon using the constant- P rule. The number of nodes for the simulated 500 random 
homogeneous point patterns has been fixed to N = 1600 and the path-loss exponent 
has been set to a = 2. Different curves correspond to k = 1 (solid), 2 (dashed), 3 
(dash-dotted); for comparison the average relative giant component is shown as the 
dotted curve. 
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Fig. 5. Probability for one-connectivity as a function of average transmission power 
upon using the local minimum-link-degree rule (curves with symbols) and the ar- 
tificial constant- P rule (curves without symbols). The number of nodes has been 
fixed to N = 1600 and the path-loss exponent has been set to a = 2. Differ- 
ent line types of the curves correspond to different random spatial point patterns: 
homogeneous (solid), multifractal (dashed) with parameters (3 = 0.4 and j = 5, 
Manhattan (dash-dotted) with parameters N x = N y = 7; a sample of 500 geomet- 
ric graphs has been used for each case. From the left to the right the symbols on 
each minimum-link-degree-rule curve stand for ngb m \ n = 3—6 (homogeneous), 3—7 
(multifractal) and 5—10 (Manhattan). 
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Fig. 6. Average relative giant component as a function of the average transmis- 
sion power upon using the iid-P rule. The number of nodes for each of the used 
500 random homogeneous point patterns has been fixed to JV = 1600 and the 
path-loss exponent has been set to a = 2. Different curves correspond to different 
parameter choices: /3i,/?2 = 0.0,0.0 (solid, constant-P rule), 0.1,0.1 (dashed), 0.3,0.1 
(dash-dotted), 0.3,0.3 (dotted). 
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(a) -LOCAL- MINIMUM-LINK- DEGREE- RULE_AD-HOC-NODE'S-VIEW() 

1 Pi ^ -fmin 

2 Mi <r- 

3 receive-mode() 

4 discovery-mode() 

5 receive-mode() 
discovery-mode() 



-i pdisc , p 

- 1 - ^ ~ i min 

2 7v; disc <- 

3 while P disc < P max and |A^ disc | < n<?6 min 

4 do Pf sc <h- Pf sc • (1 + A) 

5 SEND-HELLO(i,Pf sc ,7V; disc ) 

6 receive-hello-replys() 

7 A/f sc ^7V; disc U{ji,...} 

8 P <^max(P disc ,p) 

9 if P disc < P 

10 then SEND-HELLO(i, Pi, Mf isC ) 

11 (ki,i), RECEIVE- HELLO-REPLYS() 

12 M ^7V; disc U{fci,...} 
receive-mode() 

1 (j,Pj,Mj) <r- RECEIVE-HELLO() 

2 if i £ TV} 

3 then if \Mj \ < ngb m - ln 

4 then P, = max (Pj, P,) 

5 send-hello-reply(«, j) 

6 if |jVj| > ngb min 

7 then send-hello(«, Pi, Mi) 

8 (fei , i), . . . <r- RECEIVE-HELLO-REPLYS () 

9 Mi <-7ViU {fci,...} 



(b)-LOCAL-MINIMUM-LINK-DEGREE-RULE_SIMULATOR'S-VIEW() 

1 for i 1 to JV 

2 do P <- P min 

3 for 2 i — 1 to N 

4 do j(k) <— sort index 1 < j < N, i ^ j in increasing order of r,j 

5 P «- max (Pj, (rij(„ 9 6 min )) a ) 

6 for fc <— 1 to ngb m - ln 

7 do Pj <- max (P,-, (r ii(A )) a ) 

Fig. 7. Algorithmic implementations of the local minimum-link-degree rule: (a) the 
ad hoc node's view, and (b) the simulator's view. 
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Fig. 8. Average relative giant component as a function of average transmission 
power upon using the local minimum-link-degree rule (curves with symbols) and 
the artificial constant- P rule (curves without symbols, identical to curves of Fig. 3). 
The number of nodes has been fixed to N = 1600 and the path-loss exponent has 
been set to a = 2. Different line types of the curves correspond to different random 
spatial point patterns: homogeneous (solid), multifractal (dashed) with parameters 
[3 = 0.4 and j = 5, Manhattan (dash-dotted) with parameters N x = N y = 7; a 
sample of 500 geometric graphs has been used for each case. From the left to the 
right the symbols on each minimum-link-degree-rule curve stand for ngb m \ n = 2—7 
(homogeneous, multifractal) and 2—9 (Manhattan). 
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Fig. 9. Node degree and transmission power distributions obtained from the mini- 
mum link-degree (md.d.) and the constant-P rule. Details about the fitted distribu- 
tions are given in the main text. The number of nodes has been fixed to N = 1600 
and the path-loss exponent has been set to a = 2. Different random spatial point 
patterns are used: homogeneous (a), multifractal (b) with parameters [3 = 0.4 and 
j = 5, and Manhattan (c) with parameters N x = N y = 7; a sample of 500 geometric 
graphs has been used for each case. 



30 



